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We develop an general formalism of single enzyme kinetics in two dimension where substrates dif- 
fuse stochastically on a square lattice in presence of disorder. The dynamics of the model could be 
decoupled effectively to two stochastic processes, (a) the substrate arrives at the enzyme site in in- 
tervals which fluctuates in time and (b) the enzymatic reaction takes place at that site stochastically. 
We argue that distribution of arrival time is a two parameter function specified by the substrate and 
the disorder densities, and that it correctly reproduce the distribution of turnover time obtained 
from Monte-Carlo simulations of single enzyme kinetics in two dimension, both in absence and 
presence of disorder. The decoupled dynamics model is simple to implement and generic enough to 
describe both normal and anomalous diffusion of substrates. It also suggests that the diffusion of 
substrates in the single enzyme systems could explain the different distributions of turnover time 
observed in recent experiments. 
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I. INTRODUCTION 

Biology of life solely depend on complex chemical re- 
actions brought about by certain specific enzymes. En- 
zymes are known to facilitate reaction rates without tak- 
ing part in reactions directly [if. As catalysts, enzymes 
are receiving increasing importance from other branches 
of science. The mechanism of enzymatic activity is one 
of the most fascinating field of physical chemistry. In 
general in a chemical reactions like (A + B — > P) the 
rate of reaction depends linearly on the initial concentra- 
tion of the reactants. But in several enzymatic reactions, 
where enzyme (or catalyst) E and substrate (or reac- 
tant) S form a product P, the rate of reaction is known 
to saturate Q with increase of substrate concentration 
[S] (although the linearity remain valid for low concen- 
tration [S]), which was explained by Leonor Michaelis 
and Maud Menten 0] in 1913. In their formulation, it 
was postulated that enzyme (or catalyst) E and substrate 
(or reactant) S form a product P through a intermediate 
enzyme-substrate complex ES as 



E + S^ES^P + P 

k 2 



(1) 



Using law of mass action it was shown that the reaction 



proceeds with velocity 



V[S] 



where V and K are con- 



(IS] + K) 

stants which depend on the rates of the reaction. Over 
years this formulation, popularly known as Michaelis- 
Menten (MM) kinetics, has found applications in differ- 
ent enzymatic reactions and in several other systems 
[3, Q ■ Although Michaelis-Menten kinetic model is based 
on law of mass action its extension for ordinary or frac- 
tal spacial dimesions Q with or without disorder have 
been attempted . To incorporate effect of diffusion and 
noise these multi-enzyme systems have been studied in 
two dimension using Monte-Carlo simulation Q. 
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Recent development of experimental techniques allow 
one to do single molecule experiments [H Q , where sub- 
strates react stochastically with a single enzyme held at 
a fixed position. Thus, formation of the products here 
form a discrete time series, with interspike intervals r 
distributed as f(r). It has been shown that the rate of 
reaction (t) _1 follows MM law independent of the func- 
tional form of /(t). Thus the distribution of turnover 
time /(r), rather than MM law, could be taken as the 
characteristic property of any particular single enzymatic 
reaction system. To understand the variations of f(r) ob- 
served in experiments several theoretical models [Io|, IH| 
have been proposed. In Ref. [lOj], density of S, P, E and 
ES are replaced by their corresponding probabilities to 
explain the form of /(r). However effect of diffusion and 
noise, which are known to modify the theoretical descrip- 
tion in other low dimensional systems 12] have not been 
incorporated in these studies. 

In this article we study the single enzyme system in 
two dimension using Monte-Carlo simulation and calcu- 
late /(r). The substrate molecules S, here, are allowed 
to diffuse on the lattice and react stochastically with the 
enzyme E fixed at a particular site. We also propose an 
alternative simple method which illustrate that the dy- 
namics of this two dimensional system which explicitly 
takes care of the diffusion of the substrates can be decou- 
pled into two separate dynamics; (a) the substrates arrive 
at the enzyme site in intervals of £ which varies stochas- 
tically and its distribution is P(C)j and (b) the enzymatic 
reaction takes place at the enzyme site. We further 
argue that P(C) is a universal two parameter function 
independent of whether the system has quenched or an- 
nealed disorder, or if the diffusion of substrates is normal 
or anomalous. Given a P(C), the decoupled dynamics 
can be implemented using a single-site stochastic simula- 
tion, which greatly simplify the study by overlooking the 
full dimensionality of the system. 

The article is organized as follows. In section II we 
describe the Monte-Carlo simulation of single enzyme ki- 
netics in two dimension and calculate the distribution 
of turnover time /(r). Description of the decoupled dy- 
namics, characterization of distribution of substrate ar- 



2 



rival time P(C), and comparison of /(r) obtained from 
the decoupled dynamics model with that obtained from 
Monte-Carlo simulation are discussed in section III. In 
section IV, we show that the decoupled dynamics model 
correctly reproduce /(r) both in absence and presence 
of disorder. Finally, in scctionV we conclude the study 
with some discussions. 



II. MONTE-CARLO SIMULATION OF SINGLE 
ENZYME KINETICS IN TWO DIMENSION 

To model the single enzyme kinetics in spacial dimen- 
sion d = 2 (SEK2d) we take a (L x L) square lattice 
with periodic boundary condition where sites are labeled 
by two positive integers — 1,2, ...L. The enzyme 

is placed at the site i = L/2 = j, which is denoted as 
enzyme-site. At very other site, we have a site variable 
Sij = 0, 1, 2 which corresponds to the fact that a site is 
either vacant, occupied by a substrate S or by a prod- 
uct P respectively. The restriction that each site can be 
occupied by at most one particle captures the hardcore 
interactions present among the molecules. The density 
of the substrates [S] is denoted by p — ^f- where N$ is 
the number of substrates. 

The dynamics of the system can be described as fol- 
lows. Far from the enzyme-site the substrates and the 
products simply diffuse by making an unbiased move to 
one of the neighbouring vacant site. These moves are not 
allowed if all the neighbours are occupied. The enzyme, 
which is pinned to the enzyme-site, can either be in state 
E or makes a enzyme substrate complex ES. The en- 
zymatic reaction takes place as follows. If the enzyme 
is in state E and one of the neighbour is occupied by 
S, then E + S — ► ES by making the neighbouring site 
vacant. This process occurs with probability k\. Other- 
wise when ES is present at the enzyme-site and one of 
the neighbour is vacant, then either reaction ES — > E+S 
or ES — > E + P occurs with rates ki and k% respectively. 
The substrate S or the product P which is yield in these 
reactions goes to the vacant site. Note that, the system 
is updated in discrete time, hence the probability k that 
an reaction occurs is same as rdt, where r is the rate of 
that chemical reaction and dt = 1/L 2 . 

To calculate the rate of reaction in a multi-enzyme sys- 
tem one counts the number of products Np in time and 
then — ~^37~ gi yes the rate of reaction. Depen- 
dence of the rate on initial substrate concentration p is 
described by the MM law in a regime p <C 1. In the 
single enzyme reaction, however, product forms in dis- 
crete intervals r whose distribution is defined as /(r). 
To calculate /(r) we have done Monte-Carlo simulations 
of SEK2d and obtain the time interval r between forma- 
tion of any two consecutive products, t is measured in 
units of Monte-Carlo sweep (MCS). As shown in Fig. Q] 
the distribution of r is found to be a T-function, 
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FIG. 1: The distribution of turnover time /(t) calculated for 
p = 0.1 is a r- function [2] with (3 = 0.7, r = 17.5. Similar fit 
could be done for different p by taking (3 and to as param- 
eters. The inset shows variation of f3 and to with p. Other 
parameters are ki — 1, k2 = 0.001, fc3 = 0.03, and L = 100. 




/(r) = At 13 exp(-T/T ), 



(2) 



FIG. 2: The rate (r) 1 obey Michaelis Menten law. Inset 
shows that the inverse rate depends linearly on p _1 for low 
densities. 



where both the exponent (3 and the cut-off scale To de- 
pend on density of the substrates p. The insets of Fig. Q] 
shows these variations. Note that the distribution /(r), 
shown in Fig. [TJ differs substantially from the experi- 
mental studies, mainly because, the hardcore interaction 
here allow only four substrates in the neighbourhood of 
the enzyme which decreases the turnover probability /(t) 
for small r. Another possible reason is that the sensitiv- 
ity of the detector is insufficient in experiments to resolve 
small turnover time and the observed /(r) could be the 
result of an integrated effect. 

Since in single enzyme systems the average time inter- 
val between formation of two consecutive products is (r), 
the rate of reaction is (r) _1 . We have calculated (t) _1 
for different p and plotted it in Fig. [5] According to MM 
law inverse rate of turn over depends linearly on the in- 
verse rate of substrate density, which is described in inset 
of Fig. [21 The linear relation holds good for low densities 
as expected. Thus, it is evident that this single enzyme 
system follow MM law independent of the functional form 
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of /(t). This has been seen in the experiments of single 
enzyme systems @ recently. 

However, since the distribution of r contains more in- 
formation than the rate of the reaction (t) , in this article 
we give prime importance to /(t) and try to achieve a 
unified formulation which provide a simple method of 
calculating /(r), without doing a two dimensional simu- 
lation. This is done in the next section. 



III. THE DECOUPLED DYNAMICS AND ITS 
STOCHASTIC MODELING 

In the single enzyme systems two different dynamics 
act separately First, that the substrates diffuse (with 
hard core restrictions) through the two dimensional lat- 
tice and arrive stochastically at the enzyme site. Second, 
that the enzymatic reaction (pj takes place there. Based 
on these decoupled dynamics we introduce a stochastic 
model of single enzyme kinetics as follows. 

In this model, the substrate arrives stochastically at 
the enzyme site in intervals of £ and then the enzymatic 
reaction takes place there according to ((T|). For a given 
distribution -P(C), one can simulate this single site dy- 
namics as follows. The enzyme site has two possible 
states, E and ES. If it is in state E it continues to be in 
the same sate until a substrate S arrives there, and after 
that it is converted to ES instantly. Otherwise if the site 
is in state ES it is converted to E stochastically either 
by forming a product P with rate k% or by breaking the 
complex ES — > E+S with rate ki- However, in the mean 
time if the substrate arrives while the site is in state ES 
it falls off. Figure [3] describes this schematically. 
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FIG. 3: Schematic representation of decoupled dynamics 
model. 

Thus, given P(Q one can calculate /(r) numerically 
using the above stochastic model by noting down the 
time difference between formation of any two consecutive 
products. For the simplest model system, where sub- 
strates are always available at the enzyme site, f(r) can 
be calculated analytically. In this case, P(C) = <5(C — A), 
where A is the unit of time in discrete-time-update (here 
A = 1 MCS). Thus, the enzyme site is always in the 
state ES and form products stochastically with with rate 
fc 3 . Of course, ES breaks up with rate ki to from E + S, 
but then instantly E — > ES as other substrates are avail- 
able at the enzyme site. As expected, this results in an 
exponential distribution f(r) — A(\ — k^) T , with A be- 



ing the normalization constant. The exact form of f(r) is 
compared with that obtained from numerical simulation 
of decoupled dynamics in Fig. 0J 
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FIG. 4: Exact result (solid-line) of the turnover distribution 
f(r) = A(l - k 3 ) r for P(C) = 5(C - 1) is compared with the 
corresponding simulation of decoupled dynamics(symbols). 
Here, fci = 1, fo = 0.001, k 3 = 0.03. 

In the next section we calculate P(() for the single 
enzyme system. 



A. Distribution of arrival time P(Q 

First let us consider the system without hardcore inter- 
action among the substrates, which in turn allows each 
site of the lattice to accommodate more than one sub- 
strate. Hence each substrate in the system does a sim- 
ple random walk in two dimension. Starting from the 
enzyme site E, a simple random walk would end at a 
distance within r and r + dr in time £ with probability 

Q(C) = -±=e~ r2 /<dr. (3) 

Thus, a substrate residing in the region r and r + dr 
would reach the enzyme E with the same probability 
Q(C). Therefore the arrival probability of a substrate 
(from anywhere within the lattice) at E is 




where A depends on the upper limit of the integral and 
C is a normalization constant. Note that, for a finite 
system £ can not be larger than a typical time scale £oj 
which introduces a cut-off £o to P(£) , 

P( C )=Cerf(4)exp(-f), (5) 
vC Co 

Interaction among substrates, e.g. hardcore interac- 
tion in our study, further modify Eq. ([5]) in two ways. It 
decrease the cutoff Co which now depends strongly on the 
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density of particles p, and it also modifies the exponet of 
C which appears in the argument of erf function; 

P(C) = CerfAe X p(-f). (6) 
C Co 

In fact a depends on the exponent 7 which appears in 
the mean square displacement of a tagged particle as 

For non- interacting systems 7 = 1, resulting in a = 
7/2 = 1/2, which is consistant with Eq. (|4]). In presence 
of hardcore interaction, however, (r 2 ) gets modified by 
a logarithmic term [l3T ] in the asymptotic regime. Since 
in the enzymatic systems we are rather interested in the 
small time limit, the logarithmic term appears as an ef- 
fective exponent 7' > 1 which results in a > 1/2. 
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FIG. 5: Arrival time distribution P(C) plotted for two dif- 
ferent densities p = 0.05,0.15 with L — 100. Fhe solid lines 
here corresponds to the best fit according to Eq. ([6j with 
a = 0.64, .7 , Co = 27.5,7.8 and C = 1.05,2 for respective 
densities. For both the densities A — 0.1 

In the following we calculate numerically. Sub- 

strates, initially placed randomly, are allowed to diffuse 
on a two dimensional square lattice (with periodic bound- 
ary condition). Arrival time of any of the substrates at a 
pre-decided fixed site is recorded. The distribution P(£) 
is calculated from the difference between any two consec- 
utive arrival times. The numerical data of P(£) obtained 
for two different density of substrates p = 0.05, 0.15 could 
be fitted to Eq. © by using a and Co as fitting parame- 
ters (see Fig. [5]). An excellent fit with the numerical data 
strongly supports the theoretical form of -P(C) described 
by ©. 



enzyme problem P(C) is generically described by ([6]), a 
function with two parameters a and Co (unless otherwise 
specified we take A = 0.1 all through the article). For 
the same system we know from the Monte-Carlo simula- 
tions of single enzyme systems that f(r) is a T- function. 
First let us check that if P(C) given by © truly produce 
the same r-function in the decoupled dynamics model. 
Figure [HI demonstrates that numerical data for /(r) ob- 
tained from the decoupled dynamics model fits well with 
a r-function. 




FIG. 6: /(t) obtained from the decoupled dynamics model is 
a r-function with (/3,r ) = (.85, 40), (0.61, 34.5) for (a, Co) = 
(0.3, 50), (0.9, 50) respectively. 

Hence from the decoupled dynamics wc infer that the 
parameters of /(r) depends on P(C)- Thus, for the com- 
plete description of the model it is reasonable enough to 
characterize (/3,to) in terms of (a, Co) which is done in 
the following. Figure [7] shows the dependence of (3 on 
a for different values of Co = 15, 25, 50, 70, 100. Clearly 
from the numerical data it emerges that 



Similarly Fig. [8] describes the dependence of To on Co 
for different values of a = 0.3,0.5,0.7,0.9. From the 
figure it it is clear that To is independent of both a and 
Co below a specific time scale C*- For C > (*, however, 
functional dependence of To on Co varies with choice of 
a. We notice that the time scale C* originates from the 
product formation rate (and do not depend on other 
rates, like k\ or /C2) as 




B. Distribution of turnover time f(r) for single 
enzyme system 

We have already seen that /(t) can be obtained nu- 
merically for any given P(C) using the decoupled dynam- 
ics model discussed earlier in this section. For the single 



Thus, to compare /(t) obtained from the two dimen- 
sional simulation of single enzyme kinetics with that of 
the decoupled dynamics model it is convenient to restrict 
Co in the regime Co < C* = ^ n ^ ac ^ without loss of 

generality one may take 

Co = C = l/fc 8 . (8) 
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FIG. 7: The figure shows f3 versus a, calculated form the 
decoupled dynamics model, for Co = 15, 25, 50, 70 and 100 
(different symbols). The relation /3 = 1 — a/2 (solid- line) is 
satisfied for all values of Co- 
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FIG. 8: The figure shows to versus Co, calculated form the 
decoupled dynamics model, for different a. It is evident that 
to is independent of a for C < £*. 



C. Comparison of decoupled dynamics model with 
Monte-Carlo simulation 

From the last section we have seen that the stochas- 
tic decoupled dynamics model generates a correct form of 
/(t) for any P(C) given by Eq. To make a correspon- 
dence with the SEK2d one needs a mapping between two 
parameters of /(r), namely (/?, r ), with (a, Co) of P((). 
From the numerical characterization of P(C) described 
in Fig. [7] and we have seen that the mapping is not 
unique. The same values of (/3, To) can be obtained from 
different (a, Co) by varying the reaction rates of stochastic 
model. We take an advantage of this situation to make 
a choice ©. Now let us verify if such choice works. 

In Fig. [5] we have compared /(t) obtained from sin- 
gle enzyme kinetics in two dimension and correspond- 
ing decoupled dynamics model for two different densities 
p = 0.1,0.21. To calculate /(r) for density p = 0.1, 
which is a T-function with /3 = 0.70 and To = 17.5, we 
simulate the decoupled dynamics for this system with 
a = .6, Co = 17.5 and &3 = 0.057 as prescribed by Eq. 
©. Similar procedure is used for p = 0.21 where (3 = 0.5 
and T = 16.1. 



0.04 




Otherwise, when Co > C* = V&3> the choice of To would 
depend explicitly on the value of a (see Fig. [5]). A com- 
patible choice, which we use here, is Co = T o- Thus, to 
obtain the return time distribution /(t) which is charac- 
terized by two parameters j3 and To we make a choice of 
P(C) with parameters 



a = 2(1-/3) ; Co = r 

and &3 = 1/t (9) 



in the stochastic simulation of decoupled dynamics, 
where the last line of ([9]) corresponds to the fact that is 
constrained by l[5]). and we have used Co — r o- I n the next 
section we will check if choice ^ in decoupled dynamics 
model reproduce the /(t) same as that of SEK2d. 



FIG. 9: /(t) obtained from decoupled dynamics model is 
compared with the same obtained from SEK2d for two differ- 
ent p = 0.1, 0.21. Other parameters are ki = 1, = 0.001. 

From Fig. [9] it is clear that /(r) obtained from SEK2d 
matches reasonably well with that obtained from cor- 
responding decoupled dynamics. A small discrepancy 
which is visible for both the densities is possibly due to 
the fact that a discretized form of -P(C) and a cut-off in 
C has been used in the numerical simulation for conve- 
nience. A minor modification of Co can compensate for 
this discrepancy, /(t) with a modified Co = 18-5 for den- 
sity p = 0.1 is presented in the inset. 

From the above analysis it is evident that the formu- 
lation of decoupled dynamics model works well for single 
enzyme systems without disorder. In the next section we 
will apply the same formalism to systems with mobile or 
immobile impurities. 
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IV. SINGLE ENZYME SYSTEMS IN 
PRESENCE OF DISORDER 

In this section we discuss the effect of disorder in 
SEK2d. Disorder can be introduced in the single en- 
zyme system by introducing No impurities in the two 
dimensional lattice. In case of annealed disorder, the im- 
purities are allowed to diffuse in the system. Thus, it is 
equivalent to a late time single enzyme system without 
disorder when products are considered as diffusing im- 
purities. Compared to systems without disorder, here, 
only the cut-off time scale t is larger. So, in this sec- 
tion we choose to work in the system of quenched dis- 
order in details. The impurities in these systems are 
immobile and do not diffuse through the system. Thus 
the substrates are restricted from visiting these impu- 
rity sites, and the effective dimension of these system are 
d < 2. Simple diffusion on these fractal lattices [H[ are 
anomalous where mean square displacement is given by 
(r 2 ) ~ t 1 with 7 < 1. Note that 7 = 1 corresponds to 
the normal diffusion. Obviously, 7 depends on the dis- 
order density po = Np/L 2 for small disorder densities. 
Beyond a critical density pu > 0.4 [l5[ however the dis- 
order sites percolate and the substrates remain confined 
in certain islands. 




FIG. 10: f(r) obtained from decoupled dynamics model (solid 
line) is compared with that of SEK2d (symbols) for disorder 
density pr> = 0.1 and density of substrates p — 0.1. Here 
&3 = 0.03 for SEK2d and fc 3 = i for decoupled dynamics 
model. Other parameters are k\ = 1, fe = 0.001. 

The anomalous diffusion is expected to change the be- 
havior of arrival time distribution P(Q. In particular, if 
we ignore hardcore interaction among the substrates, Eq. 
([5]) gets modified in presence of disorder as 

P(C) = Cerf(-A)exp(-f ). (10) 

C 7/ Co 

Again, the hardcore interaction introduces a a logarith- 
mic correction term to the mean square displacement, as 
we have seen in systems without disorder in section III 
A. This correction may appear in Eq. ([5]) as an effec- 
tive exponent a ^ 7/2. Next, we ask if /(t) is still a 
T-function in presence of quenched disorder, and if the 
decoupled dynamics with a choice © reproduces correct 



/(t). It turns out that answer to both the questions 
are positive which is described in Fig. [TO] To calculate 
/(t) we did a Monte-Carlo simulation of single enzyme 
kinetics on a two dimensional square lattice (100 x 100) 
with po = 0.1 and p = 0.1. This gives rise to /(r) 
as given in © with (3 = 0.72 and r = 20. From Eq. 
© one would then expect that a = 0.57 and Co = 20 
would generate the same /(t) from decoupled dynamics 
model if fc 3 = 0.05. Calculated /(r) for both the decou- 
pled dynamics model(solid line) and SEK2d (symbols) 
are shown in Fig. [TO] Clearly /(r) for the decoupled 
dynamics model compares well with SEK2d. Note that, 
as expected, the value of a — 0.57 in case of disorder is 
smaller than a — 0.60 for the system without disorder. 

A very good agreement of decoupled dynamics with 
that of the SEK2d even in presence of quenched disor- 
der suggests that the former stochastic model provide an 
alternative simple description of single enzyme systems. 



V. CONCLUSION 

In conclusion we have studied the Monte-Carlo dy- 
namics of single enzyme system in two dimension both 
in presence and absence of quenched disorder. In both 
cases the distributions /(r) of turnover time r are found 
to be a T-function described by the exponent (3 and the 
cutoff scale tq. We argue that the dynamics of the single 
enzyme system could be decoupled to two stochastic pro- 
cesses; first that the substrates arrive at the enzyme site 
in intervals which fluctuate in time, and second that the 
reaction occurs at the enzyme site. We argue that the 
distribution of substrate arrival time -P(C) is a specific 
function ^ of two parameters a and Co- This functional 
form ([5]) is generic for systems with or without disorder. 
By choosing these parameters a and Co according to Eq. 
© , we further show that the decoupled dynamics model 
correctly reproduce /(r) obtained from the Monte-Carlo 
simulation of single enzyme kinetic models defined on a 
square lattice. 

Distribution of turnover time /(t), rather than the 
rate (r) _1 , is a characteristic feature of a specific sin- 
gle enzyme system. Following the experiments in these 
systems § several theoretical models have been proposed 
to explain the underlying cause behind the variation in 
observed /(t). In these formulations, after the enzyme- 
substrate complex breaks up by forming a product the 
enzyme returns to its normal state with a delay. Such 
a delay is considered to be essential [l(J, without which 
asymptotic decay of f(r) for large r can not be obtained. 
Further, all the reaction rates there are assumed to vary 
stochastically, which is attributed to the dynamic disor- 
der associated with several possible conformal variations 
of enzyme and enzyme-substrate complex. Alternatively, 
here we show that a simple Michaelis-Menten kinetics ([1]) 
in presence of diffusion and noise could produce the same 
/(t) observed in experiments. 
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